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Abstract. We present a microscopic analysis of the collective behaviour of the lead isotopes in the vicinity 
of 208 Pb. In this study, we rely on a coherent approach based on the Generator Coordinate Method 
including exact projection on N and Z numbers within a collective space generated by means of the 
constrained Hartree-Fock BCS method. With the same Hamiltonian used in HF+BCS calculations, we 
have performed a comprehensive study including monopole, quadrupole and octupole excitations as well 
as pairing vibrations. We find that, for the considered nuclei, the collective modes which modify the most 
the conclusions drawn from the mean-field theory are the octupole and pairing vibrations. 



PACS. 

1 Introduction 

Mean-field methods have been extensively useed to study 
a wide range of nuclear properties. There exist now several 
variants which although relying on different physical ap- 
proaches, lead to rather similar results. Among the meth- 
ods most extensively used in the last years, one can men- 
tion those based on Skyrme (zero-range) |D,p] interactions, 
or on the Gogny force (finite-range interaction) ||[|,[| , the 
density-functional approach ^ and the relativistic mean- 
held method based on an effective Lagrangian . 

One of the key ingredient in all these methods is the 
way the nucleon-nucleon interaction is treated. Either a 
density- functional or a relativistic Lagrangian is globally 
adjusted to some classes of nuclei and/or nuclear matter 
properties. Henceforth, it should be used without read- 
justment for any specific problem relevant to this class of 
properties. Mean-field methods attempt to reproduce nu- 
clear many body properties within a description which 
takes into account the one body density only. Expect- 
edly, the data accessible to these studies is mostly lim- 
ited to that corresponding to one body operators, such 
as the shape deformations or the radii. To that already 
large set one can add the total binding energy for which 
the Hartree-Fock (HF) or Hartree-Fock Bogoliubov (HFB) 
minimization principle is explicitly designed. In practice, 
other observables are sometimes also included in the fit 
of the functionals or lagrangians. A typical example is 
the effective nucleon mass which is related to the den- 
sity of single particle levels at the Fermi surface. To ob- 
tain a density similar to that deduced from the results 
of pickup and stripping reactions, the effective mass of 



the interaction should be close to the free nucleon mass. 
On the other hand, it has been shown that RPA correla- 
tions which are the simplest type of correlations beyond 
HF, increase the value of the effective mass ||.. To ac- 
count for such facts, two approaches have been utilized 
in non-relativistic studies as soon as the nuclear energy 
density approach has been shown to be flexible enough to 
reproduce nuclear properties with a good accuracy. In one 
of them, one tries to define the functional so as to avoid 
a double counting of the correlations. This implies that 
the choice of the effective parameters leaves room for the 
improvement expected from a further inclusion of RPA 
correlations. For instance, the mean-field value of the ef- 
fective nucleon mass is adjusted to a value smaller than 
the bare mass so that RPA correlations may enhance it to 
the experimental value. In the past, this has been the most 
standard approach to define effective interactions used in 
the calculation of the energy functional. Another strategy 
consists in the construction of functionals which are in a 
sense "doubly" effective since they are expected to take 
into account the effect of all correlations beyond mean- 
field on a specified set of nuclear properties. Then, for 
such properties, the mean-field calculation should yield 
results directly comparable to data. Two successful exam- 
ples of the latter strategy are the interaction of Tondeur 
et al. J10[ which reproduces most known nuclear masses 
with high accuracy and that of Brown pd| ] which yields 
single particle level schemes to be used in subsequent shell 
model calculations. 

In this work, we consider a functional belonging to 
the the first family in order to investigate how collective 
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modes modify the mean-field properties. For this analysis, 
we rely on the generator coordinate method (GCM) taking 
into account several collective variables. In the past, such 
an approach has been shown to be suited for a description 
of low lying collective states related to the quadrupole jl2| 
or the octupole [jl3| degree of freedom. A study of the cou- 
pling of the dipole and the octupole modes in 152 Sm and 
in the superdeformed band of 190 Hg (l4) has also demon- 
strated the ability of the GCM to describe in a single 
calculation the properties of the ground state, of the low 
lying collective states and of the giant dipole resonance. 
This method seems therefore to be an appropriate tool to 
study the effect of correlations beyond a mean-field ap- 
proach on the ground state properties of a nucleus. 

As test case, we have selected five Pb isotopes around 
the magic number N = 126. These nuclei are in principle 
well described by mean-field theories. The evolution of the 
two neutron separation energy at a magic number is rep- 
resentative of the shell gap. We will therefore determine 
how much it is affected by different types of collective cor- 
relations. Our aim is to be as systematic as possible within 
the GCM and to investigate the effect of all collective nu- 
clear shape degrees of freedom which are believed to be 
important for these nuclei : the monopole, quadrupole and 
octupole modes. 

Moreover, our study also includes the pairing vibra- 
tion mode. Indeed, as part of this work, we analyze the 
validity of a commonly used treatment of pairing corre- 
lations beyond the BCS theory. Since BCS correlations 
in 208 Pb vanish, one often attempts to describe them by 
means of an approximate variation after projection on the 
nuclcon numbers. In the following, we test the validity of 
the Lipkin Nogami (LN) prescription |l5|,|l6],[l7| by per- 
forming an exact projection of the BCS wave functions. A 
further study of pairing vibrations also naturally provides 
an alternative approximate variation after projection and 
allows to test the quality of the LN prescription. 

The organization of this work is as follows. In Sect.||, 
wc present the dynamical ingredients (Hamiltonian and 
collective spaces) and recall the method while introducing 
our notations. The Sect .[| discusses our results. In the last 
section, we summarize the major conclusions and indicate 
possible extensions of the present work and promising out- 
looks. 



2 Method 

In this section, we briefly review the now well established 
non-relativistic microscopic framework which leads to the 
results on lead isotopes presented in Sect. ||. It provides a 
natural setup for a treatment of i) the mean-field proper- 
ties by means of standard HF+BCS or HFB methods with 
effective Hamiltonians and ii) the collective dynamics of 
any amplitude by means of the GCM. 



2.1 The Effective Hamiltonian and the Mean-Field 
Equations 

Over the last quarter of century, it has been acknowledged 
that the separation of the nuclear binding energy at the 
mean-field level in Hartree-Fock and pairing components 
in the spirit of the energy density approach provides a 
fruitful starting point for an investigation of the static 
and collective properties of nuclei. Without requiring com- 
pletely consistency, since the Hartree-Fock and pairing 
fields remain coupled through the one-body density ma- 
trix, this approach deviates from the strict HFB theory as 
it can be performed with the Gogny || or SkP |i| interac- 
tions. Still, in nuclear physics, the energy density approach 
is often formulated in a way which preserves contact with 
the notion of effective force. Then, as it will be the case 
here, two interactions are used; one in the particle-hole 
(Hartree-Fock) channel and one in the particle-particle 
(pairing) channel. The introduction of these two effective 
forces yields additional flexibility in the phenomenological 
description of nuclear properties. 

Historically, the structure of the particle-hole interac- 
tion, which in heavy nuclei accounts for more than 99% 
of the binding energy, has been investigated first. In the 
present work, we use one of the latest versions of the well 
established line of Skyrme interactions: the SLy4 parametriza- 
tion Q designed to describe properties of nuclei (light 
and heavy), near the stability line or in the neutron rich 
region up to the drip line with a similar accuracy. Here, 
we are concerned with stable and long lived lead isotopes. 
The quality of SLy4 for such non exotic nuclei has been 
tested on numerous examples ]l9| , p0| .p|j2l| . The force which 
acts in the pairing channel and defines the corresponding 
terms in the energy density is taken from Refs. f22| , po|| . 
It is a zero-range interaction with a dependence on the 
nucleon space density. This dependence ensures that the 
particle-particle two-body matrix elements arc sensitive to 
the behavior of the single-particle orbitals mostly located 
at the surface of the nucleus. This interaction has been 
tested on physical issues sensitive to pairing properties, 
such as the evolution of moment of inertia in extended 
rotational bands as have been observed for instance in su- 
perdeformed Q and superheavy (23) nuclei. 

Along a series of isotopes or isotones, as either N or Z 
approaches a magic value, the mean-field method predicts 
a sudden collapse of pairing. It is now established that 
the magnitude of this transition overemphasizes the low- 
ering of pairing correlations which must be expected when 
the single particle density at the Fermi surface decreases 
significantly. Indeed, even in magic nuclei when the mean 
value of the order parameter vanishes, it is expected that 
fluctuations of this order parameter (in phase and magni- 
tude) should still influence nuclear properties. The impor- 
tance of this so-called dynamical pairing (as opposed to 
static pairing in situations when the HFB and HF solu- 
tions differ from each other) is one of the questions that we 
investigate below within the GCM formalism. Lipkin and 
Nogami have proposed a prescription which attempts to 
correct the mean-field for the effects of the phase fluctua- 
tions of the order parameter [|l5Ul6| , [l7|] . This still not well 
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Vibration 


Operator 


Monopole 


Qoo oc J2i r i 


Quadrupolc 


Qlm OC J2i r i Y 2m(9i,(Pi) , 771 = 0, ±2 


Octupole 





gaps parameters A n and A p : 



H P (A n , A p ) 



Table 1. Constraining operators used to construct the collec- 
tive spaces associated with the lowest isoscalar nuclear shape 
vibrations. The quantities n, Q% and tpi refer to the spheri- 
cal coordinates of the nucleon i in an intrinsic reference frame 
whose origin is at the nucleus center of mass. The Yi m are the 
standard spherical harmonics. 



understood prescription has nevertheless been shown to 
lead to global improvement (over the mean-field theory) 
for the description of nuclear properties related to pair- 
ing. We will therefore also present results obtained with 
the Lipkin Nogami prescription which is also sometimes 
presented as an approximation to the variation after pro- 
jection method (VAP). 

To conclude this section let us mention a technical 
point: when pairing correlations are present, whether at 
the static level (HFB) or because the Lipkin Nogami pre- 
scription has been implemented, we have taken care that 
no orbital of the continuum is occupied in the HFB wave 
function. In practice, this is achieved by means of a cut- 
off which excludes all orbitals whose single particle energy 
is 5MeV above the Fermi energy, from the active pairing 
space. 



2.2 The Construction of the Collective Spaces 

In this work, the set of collective wave functions {|^(9)}} 
are determined by self consistent methods (q denotes the 
value of the collective variable of interest). For the modes 
associated with isoscalar shape vibrations (monopole, qua- 
drupolc, octupole), these spaces are obtained by constrained 
HF+BCS calculations with the appropriate one-body op- 
erator as listed in Table 1 and q is the expectation value 
of the multipole moment operator q = (<P(q)\Qj m \<P(q)) 

For the vibrations in the pairing space the collective 
coordinate must be related to some gl oba l (complex) pair- 
ing gap. Since, as is discussed in Sect. ^3, the fluctuations 
associated with the phase of the gap are taken into ac- 
count by an exact projection after the HF+BCS variation, 
there only remains to describe the fluctuations of the gap 
magnitude. To build the associated collective space, there 
is no unique prescription For instance, one may consider 
generating the space by means of HFB calculations with 
a constraint on the fluctuations of the neutron and proton 



numbers (i.e. the operators AN* , fi — n,p) J2J|. 

In order to escape the technical difficulties associated 
with the two-body nature of this operator, we have adopted 
a second method. We construct the space by means of non 
constrained HF+BCS calculations using the same func- 
tional (with the SLy4 parametrization) for the mean-field 
part and replacing the pairing functional by that associ- 
ated with an auxiliary Hamiltonian depending on two real 



A 



N 



a, - 4- 



£ V a f a j -+ cc 



(1) 



where aj.„ is the creation operator of the nucleon of isospin 
t in the individual state \i and where ft denotes the time 
reversed orbital. The variation of the corresponding total 
functionals generates a set of BCS states {\<P(A n , A p ))}. 
Once this collective space has been determined, all fur- 
ther calculations, (expectation valu es and GCM matrix 
element calculations (see Sect. 2.3) ) ar e performed with 
the Hamiltonian described in Sect. |2.1| . This method has 
already been used in Ref. |2j| to investigate the influence 
of dynamical pairing on tunneling probabilities between 
the superdeformed and normal wells in 192 Hg. 

The GCM method involves a non orthonormal collec- 
tive basis. The two prescriptions for constructing the pair- 
ing collective space sketched above give should therefore 
lead to similar results (collective energy and wave func- 
tions) once the bases are large enough. In addition, the 
GCM method does not explicitly depend on which col- 
lective coordinates have been selected to span the total 
collective space. As long as the mapping is one to one, 
rather than by (A n ,A p ), we may as well choose to plot 
the GCM collective properties in terms of the coordinates 
((AN%), (ANp)) where (A) denotes the expectation value 



of the operator A in the state ^{A,- 
set of coordinates) . 



A p )) (or by any other 



2.3 The Projection and the Collective Dynamics 

The Lipkin-Nogami prescription determines a BCS wave 
function which, in priciple, takes into account the pairing 
fluctations associated with the phase 4> of the pairing. This 
wave function can be used directly to evaluate the expec- 
tation value of any operator as proposed in the original 
articles of Lipkin and Nogami. In this work, such values 
will be refered to as L.N. One can also consider that this 
BCS state is an approximation of the exact VAP state. 
In order to test the validity of this assumption, one must 
then extract from the LN-BCS state its component with 
correct (i.e. No) particle number by means of a projec- 
tion Pn - Henceforth, the results associated with such a 
two-step approach will be labeled L.N. (proj.) 
The projection 



(2) 



involves an integral over <f> with the weight factor e~^ N ° . 
From the mean-field wave functions { \$(q)) } associated 
with the value q of the collective variable, we build the 
states 



(3) 
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which form the non orthogonal projected collective basis 
{|^jv (<?))} The expectation value £(q) of the Hamiltonian 
in this basis is the projected deformation energy curve: 



($ n m\h\*nM) _ mq)\Hp No \m) 



(*jv D (?)l*JVo(g)> 



($(q)\P No \$(q)) 



(4) 



Several such curves are discussed in Sect|. The value at 
the minima of £(q) yields an energy corresponding to a re- 
stricted VAP (a variation limited to the subspace spanned 
by the collective variable q). The quality of this upper 
bound to the VAP energy depends on the relevance of the 
collective space for dynamical pairing correlations. 

In a next step, we consider a more general N-body 
wave function defined as a linear superposition of pro- 
jected HFB states with an unknown weight function fj(q): 



(5) 



In this definition, the label j recalls that several states 
are obtained corresponding to the correlated ground 
state and to the collective excited states. Within the GCM 
method, the function fj is determined by a variation of the 
total energy Ej 



(6) 



with respect to the function fj(q)- This leads to the Hill- 
Wheeler equations p^] . When particle number projection 
is imposed on each basis state, the kernels X, M and H of 
the integral operators entering this equation are given by 



H(q',q) 



-i<f>N 



<*(«')! 




The evaluation of these (q',q) dependent kernels is the 
most time consuming numerical part (in fact, there are 
four of them since neutron and proton numbers have to 
be conserved separately). They involve a double integral 
(one for N and one for Z) of matrix elements of one-body 
and two-body operators between all possible states of the 
collective basis {^(q))}- The formulae useful for the calcu- 
lation of the sets of matrix elements are given in Ref. |^7j . 
Other kernels involving multipole moment operators asso- 
ciated with the various deformations are also needed for 
the computation of the shape properties of the collective 
GCM states. The reference just quoted also recalls how 
the solution of the Hill- Wheeler equation can be reduced 
to the diagonalization of an hermitian matrix computed 
from the kernels X, AT and Ti. 

The functions /j's do not form an orthogonal set. Using 
the overlap kernel, one can define the functions g^s 



g j (q) = (XV 2 .f j )(q) 



(8) 



which are orthonormal and can be interpreted as collec- 
tive functions in the usual sense. In Eq. (|8|), the notation 



X 1 ! 2 stands for the integral operator whose square is equal 
to the integral operator with kernel X(q' 1 q). In the next 
section, we present some functions gj associated with the 
quadrupole and octupole vibrations. 

In the case of quadrupole vibrations, the GCM is solved 
separately in each representation of the permutatio n group 
of the three intrinsic axes. As discussed in Ref. pq , p9[ , 
the two independent representations of this group have 
different angular-momentum-parity (J*) contents: in the 
fully symmetric one, the major component of the collective 
wave function corresponds to + spin and parity, while, in 
the two dimensional representation, the collective N-body 
wave function is predominantly of the 2 + type. Very re- 
cently, we have contructed an exact angular momentum 
projection code [j30t and we have quantitatively checked 
the quality of this approximate spin projection. We have 
found that it is accurate whenever the extension of the 
collective function g does not extend beyond the value 
(Qw) = 4b of the quadrupole deformation. Since the in- 
fluence of the spherical magic proton number Z — 82, in- 
sures that g is non zero only in the close vicinity of the 
(Q20) = point, this is the case for the g functions of the 
lead isotopes. 

Although our HFB and GCM codes allow triaxial de- 
formations, the results presented below are based on col- 
lective bases including only prolate and oblate states. This 
is justified by the small extension of the collective wave 
functions. Nonetheless, all the GCM calculations presented 
below mix states located on the six semi-axes correspond- 
ing to the triaxiality angles 7 = 0, ±60°, ±120° and 180°. 
We have tested the quality of this approximation by a 
performing full triaxial calculations on a few cases and we 
have checked that our conclusions do not depend on this 
degree of freedom. 

The mean-field (and Lipkin-Nogami) wave functions 
spanning the collective spaces are constructed from an en- 
ergy density functional whose particle-hole and particle- 
particle components are defined by means of two distinct 
effective density dependent two-body forces. For the sake 
of simplicity, we have presented the standard GCM for- 
malism with a single effective Hamiltonian. In pract ice, 
we have adapted the GCM to the choice made in Sect. 2.1 
and defined a generalized energy density to be used in 
the computation of the energy kernel TC(q',q) given by 
Eq. (Q) . This involves both a prescription for handling the 
density dependence of the force and an extension of the 
energy density formalism to the non-diagonal matrix ele- 
ments which appear in the integral in Eq. ([?]) . We use the 
method described in Ref. |2£| (for an alternative choice, 
see Ref. |n) for instance) . The formulae for the two terms 
of the extended energy density are given in Ref. Q and 
the implications and potential problems associated with 
this formulation are discussed in Refs. |32|,E3. 



3 Results 

We consider successively, the influence of correlations on 
the ground state properties and the collective excitations 
of five lead isotopes. 
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Fig. 1. Differences between the calculated total binding en- 
ergy of the ground state of the five lead isotopes and a BCS 
reference. The labels Qm, i = 0, 2, 3 refer to GCM calculations 
taking into account the different types of shape collective exci- 
tations reported in Table [l], while "Pair." stands for the pairing 
vibrations. Except for the curves BCS and L.N., all results in- 
volve an exact projection on N and Z. 



3.1 Ground State Correlation Energies 

In Fig. 0, for each of the isotopes 204 Pb, 206 Pb, 208 Pb, 
210 Pb and 212 Pb, we have plotted the difference between 
the binding energies of the ground state calculated by 
GCM on several collective spaces and a BCS reference 
value. In practice, for the five nuclei, the proton part of 
this reference wave function reduces to a Slater (HF) de- 
terminant. This is also true for the neutron component in 
208 Pb. In addition, in Fig. |], we have plotted the difference 
between the same reference and the binding energies cal- 
culated by means of the Lipkin Nogami method with and 
without a projection on the neutron and proton numbers. 

The ordinate scale of the figure shows that, irrespective 
of the nature of the collective vibration, the correction to 
the total binding energy is small (about 3 MeV). Consid- 
ering the effective character of the two-body interactions 
defining the energy density, this is a satisfactory result. 
As we discussed in the introduction, one can consider ei- 
ther that, to first order effects of correlations can be ef- 
fectively included into the parametrization of the force, or 
conversely that the presently available forces have to be 
readjusted so as to underbind slightly with respect to data. 
This would leave room for the small corrections associated 
with correlations to be computed on top of the mean-field 
solution, as we have done. In Ref. much larger values 
of correlation energies have been found. However, as was 
pointed out in this reference, the method based on a sum 
rule approximation to the RPA suffers from an ambiguity 
related to the double counting in RPA formulas for the 
correlation energies. The GCM method is free from such 
a problem. 



The relative magnitudes of the correlations tell about 
the nuclear sensitivity to various collective modes. For lead 
isotopes the octupole mode is the most efficient among the 
three nuclear shape vibrations analyzed in this work. One 
also notes that even for nuclei close to N and Z magic shell 
closures, pairing vibrations generate the largest correction. 

The L.N. and even more the the L.N. projected curves 
show a larger effect for the doubly magic nucleus 208 Pb. 
For the latter curve, the difference with BCS increases 
from about 0.65 MeV at 204 Pb and 212 Pb to 1.4 MeV 
at 208 Pb. A possible interpretation is that the LN pre- 
scription which is designed to correct for vanishing pairing 
fluctuations in the HF+BCS solution, is more effective in 
208 Pb when the single particle energy gap at the Fermi 
surface is larger. With respect to the second curve (L.N.), 
the monopole, and octupole vibrations introduce an addi- 
tional lowering of the energy whith an almost N indepen- 
dent value equal to 0.2 MeV and 1.4 MeV respectively. 
By contrast, the curve associated with pairing vibrations 
shows an overall decrease versus N of the binding energy 
correlations : over five isotopes it is reduced by about 0.6 
MeV. 

In Fig. ||, we show the same results when the stan- 
dard LN energy is taken as a reference. This presentation 
illustrates the improvements brought about by an exact 
projection as compared to the approximate prescription 
introduced by Lipkin and Nogami. 

When the LN prescription does not yield a good ap- 
proximation of the binding energy of the particle projected 
VAP solution, one could hope it to be a practical method 
to construct a BCS state close to the VAP intrinsic state. 
Then, a projection of this BCS wave function should give 
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Fig. 2. Differences between the calculated total binding en- 
ergy of the ground state of the five lead isotopes and the Lipkin 
Nogami reference. The labels Qio, i = 0, 2, 3 refer to GCM cal- 
culations taking into account various types of shape collective 
excitation as described in Table |l[ while "Pair." stands for the 
pairing vibrations. Except for the curves BCS and L.N., all 
results involve an exact projection on N and Z. 
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(1.85, 1.17) 
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(1.36, 1.17) 


(1.92,1.12) 


-0.21 


210 


(1.92, 1.17) 


(2.37,1.13) 


-0.84 


212 


(2.59, 1.17) 


(2.58,1.50) 


-0.28 



almost the same energy as a VAP calculation. Our re- 
sults show that the validity of this property of the LN 
prescription is at least nucleus dependent. The minimum 
of the energy surface obtained with the projected BCS 
states spanning the pairing collective space provides an 
upper bound of the VAP energy. From the results given 
in Table |[ it is seen that for the lighter lead isotopes, 
it is more than lMeV below the energy associated with 
the projected Lipkin Nogami state. The same table shows 
that the lowest binding energies correspond to projected 
BCS states with particle number fluctuations larger than 
those predicted by the LN prescription. For 208 Pb and iso- 
topes above, the differences between the projected LN and 
approximate VAP solutions are smaller. 

The discontinuity of the two- neutron separation energy 
S211 at a shell closure predicted by mean-field methods is 
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Fig. 3. Comparison between the experimental two neutron 
separation energies of the lead isotopes and the values calcu- 
lated within the BCS, projected LN and GCM methods. In the 
latter case, the collective space allows fluctuations of the gap 
parameter. 
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Table 2. Neutron and Proton particle number fluctuations 
of the BCS states associated with the Lipkin Nogami solution 
and the minimum of the variational space constructed for the 
study of pairing vibrations (see text). The quantity AE give 
the difference between the energy of the latter state and the 
LN state after particle number projection has been effected (i.e. 
L.N. proj.). 
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Fig. 4. Quadrupole deformation energy curves of 208 Pb cal- 
culated with the Lipkin-Nogami state without or with exact 
projection on the neutron and proton particle number. In the 
latter case, the curve in the subspace containing the + compo- 
nent of the deformed function is also given. The solid straight 
segments give the GCM collective energy in this latter sub- 
space, while the dotted segments give the energy in the sub- 
space containing the 2 + components of GCM the collective 
functions. 



often larger than found experimetally. This effect can be 
seen in Fig. [|. The magnitude of the drop of the BCS S^n 
value between 208 Pb and 210 Pb exceeds data by approx- 
imately 1.8 MeV. Although not shown in the figure, we 
have checked that the differences between the BCS and 
experimental S^n's are not larger than 0.2 MeV for mass 
number A smaller than 204. Therefore the larger discrep- 
ancy at A = 208 cannot be ascribed to the asymmetry of 
the effective force SLy4 and suggests a deficiency of the 
HF+BCS method for magic or near magic nuclei. This as- 
sumption is confirmed by the fact that the S^n's calculated 
with the GCM on a collective space allowing the fluctua- 
tions of the pairing gap leads to an improvement by about 
0.6 MeV (i.e. 33%). By contrast, the Lipkin-Nogami pre- 
scription does not improve the mean-field results although 
it is supposed to approximately take into account the ef- 
fect of these pairing vibrations. We have not plotted the 
curves with the energies resulting from the GCM asso- 
ciated with shape vibrations. Indeed, they are not much 
different from the LN curve. 



3.2 Collective Functions and Excitation Spectrum 

The Z = 82 and N — 126 magic numbers have a definite 
influence on the individual and collective properties of the 
five isotopes investigated in this work. They strongly favor 
a mean-field solution corresponding to a HF (no pairing) 
spherical wave function. The Lipkin Nogami prescription 
which transforms the HF states into HF+BCS ones by 
enforcing a smooth decrease of the occupation probabili- 
ties from the hole to the particle states does not modify 
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substantially the deformation properties of the collective 
energy surfaces. As it can be seen in Fig. ||, the LN energy 
curves of 208 Pb, with and without exact particle number 
projection, exhibit a well marked minimum at the spher- 
ical point. The curves are very similar, the main effect of 
the projection being a downward shift by about 0.5 MeV. 
When an approximate angular projection on the + state 
is performed, the bottom of the well is flattened over an 
extension of ±4b. Qualitatively, such an effect is expected 
on general grounds. The LN wave function at the spherical 
point is a pure + state and is unaffected by projection 
while a deformed mean-field wave function is a mixture 
of different spin components. The energy of the + com- 
ponent is lower than that of all other spins and therefore 
than that of the deformed HF+BCS wave function. 

The GCM energies of the lowest collective quadrupole 
(approximate) + and 2+ states of 208 Pb are also plotted 
in Fig. ||. The first excited state is a 2 + . A comparison 
with the GCM monopole calculation shows that the lowest 
+ collective state in the spectrum of 208 Pb is the one 
shown in Fig. [|, the energy of the first excited monopole 
GCM state being at 13.4MeV. This latter value which does 
not vary with N over the five isotope series, agrees with 
monopole giant resonance data |p4 35 . This nice result is 
related to the nuclear matter incompressibility of the SLy4 
force, 220 MeV, a value which is known to be coherent 
with experiment ||,|37|f|] . 

Fig. H displays the GCM probabilities (i.e. \g\ 2 where 
g(q) is defined by Eq. (||)) for the first + and 2+ states. 
The former is peaked at the spherical point. Its width 
(« 4b) corresponds to the quadrupole extension of the 
bottom of the angular momentum projected energy curve 
shown in Fig. f|. The collective probability distribution of 
the first 2 + GCM state displays a node at Q = 0. Its 
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Fig. 5. GCM collective functions associated with quadrupole 
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function of the lowest energy state in the subspace containing 
the + component of the total wave function, while the dashed 
curve gives the same information for the subspace with the 2 + 
components. 
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Fig. 6. Octupole energy deformation energy curves calculated 
with the Lipkin-Nogami state, without or with exact projection 
on the particle number. In the latter case, the curve in the sub- 
space containing the positive and negative parity components 
of the deformed function is also given. The solid (resp. dotted) 
segments give the GCM collective energy in the positive (resp. 
negative) parity subspace of the GCM collective functions. 



extension is almost symmetric on the prolate and oblate 
sides with a slight preference for the oblate deformation. 

Four octupole deformation related energy curves for 
208 Pb are plotted in Fig. || against positive expectation 
values of the constraining operator Q30 (see Table 0). 
Since this operator is parity negative, the negative ab- 
scissae part of each curve can be deduced by symme- 
try with respect to the ordinate axis. The LN and LN 
proj. curves display a minimum at the spherical point 
(Q30) = 0- The projection merely increases the binding 
energy by 0.6MeV. An additional projection of the wave 
function on the total positive (resp. negative) parity sep- 
arates the subspace which contains the + (resp. the 3~) 
component of the LN particle projected state. The LN 
and LN proj. curves only differ significantly for small val- 
ues of (Q30} (< 2.10 3 fm 3 ), namely as long as the overlap 
between the LN particle-projected wave function and its 
parity reversed is non zero. The positive parity curve has a 
minimum for (Q30) ~ 10 3 fm 3 . Below, we will see that this 
number provides a measure of the octupole fluctuations 
in the + GCM state. We recall that symmetry garan- 
tees that the average value of (Q30) vanishes exactly for 
all GCM states and, in addition, that the negative parity 
projected energy curve diverges at (Q30) = 0. 

The GCM energies given in Fig. |6| display an alternat- 
ing sequence of positive and negative parity states. The 
first octupole GCM excited state should correspond to 
the lowest experimental 3~. Its probability distribution 
(\g(q)\ 2 ) and that of the 0+ are plotted in Fig. 0. They 
give an information on the magnitude of the fluctuations 
in the corresponding collective state. Returning to Fig. ||, 
one checks that each of these two wave functions displays 
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Fig. 8. Evolution of the lowest excited collective states of five 
lead isotopes near the magic nucles 208 Pb. The figure compares 
data with results of GCM including an exact projection on 
particle (N and Z) numbers and an aproximate projection on 
angular momentum. 



the features expected for the quantum ground state in the 
projected Lipkin Nogami curve with the associated parity. 

Up to now, the discussion in this subsection which cen- 
tered on deformation energy curves and collective prob- 
ability distributions has been concerned only with the 
doubly magic nucleus 208 Pb. The same curves calculated 
for the four other isotopes are qualitatively similar while, 
they yield significantly different collective energy spec- 
tra. Although one can discern evolutions of the energy curves 
versus N which are at the origin of the differences found in 
the GCM spectra, they turn out to be neither spectacular 
nor linked to any specific part of the collective space. For 
this reason, we chose not to plot them; they only display 
a smooth modification of the same global pattern along 
the isotope series. 

In Fig. ||, the collective GCM spectra for the three 
first excited states (0 + and 2+ quadrupole, 3~ octupole) 
are given. Qualitatively, the overall pattern of data is well 
reproduced. The positive parity levels are pushed to high 
energy as one crosses N = 126. As a consequence, the low- 
est collective state in 208 Pb is a 3~. As one moves away 
from the shell closure, the collective spectrum returns to 
the usual midshell pattern with the 2 + state being the 
lowest. However, the collective states whether quadrupole 
(0 + and 2 + ) or octupole are found by the GCM at ener- 
gies between 1.5 to 2 larger than observed. Although the 
formalism of the GCM does not distinguish potential from 
collective inertia effects, if we were to adopt the concepts 
of the standard collective model to analyze our discrep- 
ancy with data, we could say that the results are con- 
sistent with the calculated moments of inertia of the low 
energy shape vibrations being twice too small. In the next 
section, based on the experience derived from Refs. ||@ 
and from the existence of significant dynamical pairing ef- 



fects evidenced by the results discussed in Sect. 3.1 



we 



put forward possible explanation for the GCM results. 

With the Skyrme interaction SLy4, the quadrupole col- 
lective transition matrix element 208 Pb B{E2) t =7.1 W.u. 
is close to the data (8 W.u.). The agreement is as good for 
the octupole transition since we find B{E%) f=18.7 W.u. 
as compared to the experimental value equal to 32 W.u. 
[ ^4|j4l| , ^2] , ^5| . We note however that these results are very 
sensitive to the interaction since a similar calculation with 
the SGII interaction yields 14 W.u. and 58 W.u. respec- 
tively. 



4 Conclusions and Perspectives 

We have performed an analysis of the collective correla- 
tions and modes of excitations in the vicinity of the magic 
nucleus 208 Pb. One of the characteristics of this work is 
an attempt to maximize the coherence of the description 
of the mean-field and collective effects for any amplitude. 
This is done by means of the GCM solved in collective 
bases. Depending on the type of collectivity, the GCM 
bases arc obtained from constrained calculations (shape 
vibration) or derived from a set of auxiliary Hamiltonians 
(pairing vibrations). Another element of consistency is the 
use of the same effective particle-hole (p-h) and particle- 
particle (p-p) interactions in the complete sequence of cal- 
culations. These two forces had been determined in earlier 
analyses of ground state and high spin properties by means 
of mean-field methods. 

We find that correlation effects on the total binding 
energies are small. A posteriori, this result justifies the 
perturbative-like calculation of such correlations on top 
of the mean-field solution as done in the present work. It 
also supports attempts at the definition of a force to be 
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used exclusively in mean-field calculations and designed 
to effectively incorporate the small effect of correlations 
|lO| , p]| . The latter type of force is a practical tool for the 
extensive nuclear mass tabulations required for a global 
analysis of astrophysical nucleosynthesis. 

Furthermore, we have found that pairing vibrations 
induce a softening of the rapid variation of the S2N curves 
at the crossing of the magic shell iV = 126. Such accidents 
of the nuclear mass chart have a significant influence on 
the understanding of nuclear abundance peaks resulting 
from high neutron fluxes in stellar events such as class II 
supernovae explosions. A check of the persistence of this 
effect at the points where the r-process path crosses the 
Z = 82 magic line (on the TV-rich side) appears therefore 
desirable. 

Our work also provides an illustration of the validity of 
the Lipkin Nogami prescription. The LN method should 
give good approximations of both the VAP wave function 
and energy. More restrictively, it is considered that even 
when the LN energy does not reproduce well the VAP 
energy, an exact projection performed on the LN wave 
function does so f43|| . This implies that VAP is well ap- 
proximated by a projection after the solution of the LN 
equations. Our calculations show the limitations of these 
two assumptions. We find that even after projection the 
discrepancy on the binding energy is nucleus dependent 
and can be larger than l.OMeV. Our calculation provides 
only a lower bound of this discrepancy since the varia- 
tional space used in the present work only covers a fraction 
of the space available to a full VAP calculation 

We find that, for the five lead isotopes, the GCM col- 
lective vibrations associated with our variational spaces 
only generate small fluctuations around the spherical shape 
and do not modify significantly the expectation values of 
multipole moments. The behaviour of the low energy col- 
lective spectrum is qualitatively correct. In particular the 
relative positions of the + , 2 + and 3~ states along the 
isotope sequence is well reproduced. We also find a rea- 
sonable agreement on electromagnetic transition matrix 
elements. However, the energy scale of our spectrum is 
too large by a factor between 1.5 and 2 compared to data. 
The reason for this discrepancy is not yet clear. This is 
all the more disconcerting that RPA calculations on the 
Hartree-Fock solution using the same interaction provide 
a much better agreement with the low energy collective 
spectrum of 208 Pb Q. For this nucleus, we have addi- 
tionally performed a GCM calculation of quadrupole and 
octupole vibrations using the SGII interaction for which 
RPA calculations are also available . We find a similar 
pattern of fair agreement on transition matrix elements 
and discrepancy on collective spectrum energies. 

From these results, some directions for further work 
can be outlined. First it appears useful to extend the cal- 
culation to other regions of the mass table to test the 
importance of pairing vibrations at shell closure. In par- 
ticular one should consider other proton magic numbers 
with a focus on the neutron rich side at places of interest 
for the r-process. Second, one sould extend the formal- 
ism with a coupling between the vibrational modes which 



have been detected to be especially active. For the lead 
isotopes, the first such investigation should be devoted to 
a coupled calculation of pairing and octupole vibrations. 

An other point which needs clarification is the discrep- 
ancy on the low energy collective spectrum density. The 
problem does not seem related with the interaction since 
RPA provides a rather good agreement with data. We note 
also that with the same forces, GCM and RPA results on 
giant resonances agree together (and also with data). A 
first reason for the differences on low energy collectivity 
may be due to the use of a GCM basis of wave functions 
including dynamical pairing while RPA is performed on a 
Hartree Fock solution. 

A discussion of this problem can be found in J|9|,[l0| . 
Indeed, in cranking type calculations, pairing reduces col- 
lective mass parameters and leads to a decompression of 
the spectrum. In such a case, one would have to recon- 
sider which fraction of the pairing must be assigned to 
static (i.e. BCS or HFB) or dynamic (fluctuations of the 
magnitude of pairing gap) correlations. Another source 
for the different GCM and RPA results may be due to the 
definition of the GCM collective spaces. To study shape vi- 
brations, we have used a collective basis generated by self 
consistent calculations with constraining operators pro- 
portional to multipole moments. A description of the low 
energy mode could require other radial form factors more 
peaked at the nuclear surface than the present r™ form 
factors. The RPA transition form factors may provide in- 
dications for more appropriate constraining operators. 

As a conclusion, the present work has demonstrated 
the interest of our approach and also shown some of its 
limitations in its present form. Nevertheless we believe 
that it establishes that the GCM performed on many-body 
bases including exact symmetries such as particle number, 
parity (as is done here) but also angular momentum has 
strong potential for the description of nuclear properties 
which remains largely to be explored. 
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